Measurement of dynamic Stark polarizabilities by analyzing spectral lineshapes of 

forbidden transitions 



D. R. Dounas-Frazer, 1 '^ K. Tsigutkin, 1 A. Family, 1 and D. Budker 1,2 

department of Physics, University of California at Berkeley, Berkeley, CA 94720-7300 
^Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720 

(Dated: September 30, 2010) 

We present a measurement of the dynamic scalar and tensor polarizabilities of the excited state 
|5d6s 3 Di) in atomic ytterbium. The polarizabilities were measured by analyzing the spectral line- 
shape of the 408-nm 6s 2 1 So — > 5d6s 3 Di transition driven by a standing wave of resonant light in 
the presence of static electric and magnetic fields. Due to the interaction of atoms with the standing 
wave, the lineshape has a characteristic polarizability-dependent distortion. A theoretical model was 
used to simulate the lineshape and determine a combination of the polarizabilities of the ground 
and excited states by fitting the model to experimental data. This combination was measured with 
a 13% uncertainty, only 3% of which is due to uncertainty in the simulation and fitting procedure. 
The scalar and tensor polarizabilities of the state |5d6s 3 Di) were measured for the first time by 
comparing two different combinations of polarizabilities. We show that this technique can be applied 
to similar atomic systems. 

PACS numbers: 32.90.+a, 32.70.Jz, 32.60.+i 



I. INTRODUCTION 

Static (dc) and dynamic (ac) electric dipole polariz- 
abilities determine the response of neutral particles to 
applied electric fields. They are related to a host of 
atomic and molecular quantities, including the dielec- 
tric constant, refractive index, and Stark shift [IHI], an d 
are an important consideration for many current atomic, 
molecular, and optical physics experiments [5]. For ex- 
ample, polarizabilities play a vital role in the production 
of light traps for quantum information processing appli- 
cations 6 . In the context of optical atomic clocks, Stark 
shifts constitute an important systematic effect that must 
be controlled [7H10). Similarly, Stark shifts also con- 
tribute to systematic effects in atomic parity violation 
(APV) measurements [TlT414] . Hence the determination 
of polarizabilities is a priority for high-precision atomic 
physics. 

Present experiments typically rely on theoretical calcu- 
lations of electric dipole polarizabilities [15lll7| . Several 
methods for measuring polarizabilities also exist. Early 
schemes involved the deflection of atoms in an inhomo- 
geneous electric field [3j. More recent techniques include 
absolute frequency measurements [5], atom interferome- 
try [T8 - 20 , and a technique that uses light force [5T] [55] . 
However, these methods typically provide information 
about the polarizability of an atom in its ground state [5] . 
Therefore, they are inappropriate for high-precision ex- 
periments where the polarizabilities of excited states are 
relevant. 

As part of an ongoing investigation of parity violation 
in atomic ytterbium (Yb) jT3J|T3] , a scheme for measuring 
a combination of polarizabilities of the ground and an ex- 
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cited state of Yb was developed [53] . The scheme involves 
the simulation and measurement of the spectral lineshape 
of a forbidden electric dipole transition driven by a stand- 
ing wave of light in the presence of a dc electric field. 
Due to the standing wave, the ac Stark shifts of the up- 
per and lower states introduce a polarizability-dependent 
distortion in the lineshape, a phenomenon which was first 
observed and characterized during a search for APV in 
cesium [TT]. The difference of polarizabilities of the two 
states is treated as a variable parameter in the simula- 
tion and is measured by fitting the simulated lineshape 
to experimental data. We call this scheme the Lineshape 
Simulation Method (LSM). The LSM can be generalized 
to an arbitrary atomic species. 

In this paper, we present the next generation of the 
LSM. The numerical procedures accommodate for a 
broad domain of values of input parameters, e.g., the 
intensity of the standing wave. In addition, the indepen- 
dent dimensionless parameters that determine the line- 
shape have been explicitly identified, thus facilitating er- 
ror analysis. In general, the LSM is compatible with 
a variety of atomic species and field geometries. The 
LSM is sensitive to the difference of polarizabilities of 
the atom in its ground and excited states. Neverthe- 
less, this method may yield unambiguous measurements 
of the vector and tensor components of the excited state 
polarizability, as we will show. In this sense, we present 
a versatile method for measuring the excited-state polar- 
izabilities of atoms. 

We demonstrate the LSM using the 408-nm 
6s 2 1 So — > 5d6s 3 Di transition in atomic Yb. Whereas 
the previous results [23j were obtained in the absence 
of a magnetic field, the present work uses a magnetic 
field to isolate Zeeman sublevels of the excited state. 
The ac Stark shifts of the sublevels are characterized by 
different combinations of scalar, vector, and tensor po- 
larizabilities. This approach allows for the unambiguous 
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determination of these polarizabilities. The ac scalar and 
tensor polarizabilities of the excited state |5d6s 3 Di) in 
Yb are measured for the first time. Due to improvements 
in the experimental apparatus, the signal-to-noise ratio 
of the observed lineshape is an order of magnitude larger 
than for the previous implementation. In the present 
work, the statistical error introduced by the LSM is 
negligible compared to the systematic uncertainty of the 
experiment. 

This paper is organized as follows. In Section [TTJ we 
introduce our conventions for the polarizabilities and the 
Stark shift. The theoretical model, numerical procedure, 
and results of the simulation are discussed in Section [TTTJ 



In Section IV we apply the LSM to the Yb system and 
present the results. Finally, a summary of the results and 
an outlook for future experiments are given in Section [V] 



II. ATOMIC SYSTEM 

Throughout this work, we consider a Stark-induced 
transition between two atomic states of the same par- 
ity. The transition is induced by applying a uniform 
dc electric field Ed c - We assume that the transition is 
driven by a standing wave of light formed by two counter- 
propagating waves with the same polarization traveling 
in the ±k directions. In this case, the electric field of the 
light is given by 



E(r,i) = E(r) cos(wt) e, 



where 



E(r) — E cos(k • r)e 



2/2 



(1) 



(2) 



Here Eq, k, oj, and e are the amplitude, wave-vector, 
angular frequency, and polarization of the electric field, 
respectively, is the radius of the standing wave, and 
rj_ = |r — k(k • r)| is the perpendicular distance from 
the center of the standing wave. In addition to these pa- 
rameters, we define the wave-number k = |k| = oj/c and 
the wavelength A = 27r/fc, where c is the speed of light. 
We assume that Eq > since the overall sign of the field 
can be incorporated into the polarization e. Equation |2| 
is appropriate for the case of a light field with a Gaus- 
sian profile. The discussion is limited to optical frequen- 
cies. In this regime, E(r, t) is uniform over atomic length 
scales. In order to study the influence of the magnetic 
structure of the transition, we also assume the presence 
of a uniform dc magnetic field B. The quantization axis 
(z-axis) is chosen so that B = B z for B > 0. 

The dynamics of an atom in the presence of the exter- 
nal magnetic and electric fields described above is gov- 
erned by the total Hamiltonian 



H = H + Hi (t) , 



(3) 



where H Q = H A + H z + H dc and H^t) = H &c {t) 
are the time-independent and time-dependent parts of 



H. Here H A is the atomic Hamiltonian, Hi — —fi ■ B 
is the Zeeman Hamiltonian, H dc = — d • Ed c and 
H ac (t) = — d • E(r, i) are the dc and ac Stark Hamilto- 
nians, and fj- and d are the magnetic and electric dipole 
moments of the atom, respectively. We assume that B, 
Edc, and E(r, t) are sufficiently weak that Hz, H dc , and 
i? ac can be treated as successive perturbations to H A . 

Let \l) = \jJM) and £ A (£) = £ A (^J) represent the 
degenerate eigenstates of the atomic Hamiltonian Ha 
and their corresponding energies, respectively. Here 
J is the total angular momentum quantum number, 
M E {J, J — 1, . . . , — J} is the magnetic quantum num- 
ber corresponding to the projection of the total angular 
momentum along the z-axis, and 7 is a set of other quan- 
tum numbers. Then, to lowest order in the perturbing 
fields B and E dc , the eigenstates of Hq are 



\i) = \e) + j2\n 



{t'\H dc 



with corresponding energies 



£(£) = £ A (£)+£ z (£)+£ dc (£). 



(4) 



(5) 



Here £z(£) and £ dc (£) represent the Zeeman and dc 
Stark shifts, respectively. The Zeeman shift is given 
by £%{£) — ggfioBM, where ge is the Lande factor of 
the state \£) and fj,Q is the Bohr magneton. Through- 
out this work, we assume that B is sufficiently strong 
to completely isolate the Zeeman sublevels of \£). The 
dc Stark shift is given by £ dc (£) — — (V^o^-Eyc where 
E dc = I Ed c I is the magnitude of the dc electric field, and 
a'j is the dc polarizability of the atom in state \£). To 
derive Eq. Q , we neglected mixing of atomic eigenstates 
due to the magnetic field. 

The atomic energy levels are also shifted by the ac 
Stark shift [THI], which is induced by the dynamic field 
E(r, t). We assume that the frequency of the stand- 
ing wave satisfies uj s» oj qs , where uj ag is the resonant 
frequency of the electric-dipole (El) transition from the 
perturbed electronic ground state \g) to a perturbed ex- 
cited state \a). Thus it is appropriate to make a two-level 
approximation that involves neglecting dynamic interac- 
tions between states other than \g) and \a). However, 
such an approximation can only account for ac Stark 
shifts that arise due to mixing of the states \g) and \a) 
with each other. To address this situation, we modify the 
energy of the perturbed ground state as follows: 



£{g)^£{g)+£ &c {9), 



where 



£ ac (g) 



agE(v)\ 



(6) 



(7) 



is the ac Stark shift of \g) due to mixing of \g) with 
states other than |a), and a-g is the corresponding ac 
polarizability of the state \g). An analogous modification 
is made to the energy of the perturbed state \a). 
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FIG. 1: Energy-level diagram. Shown are the energy eigen- 
states of an atom and the electronic transitions relevant to 
the LSM. 



We further assume that the unperturbed ground and 
excited states \g) = \j g J g M g ) and |o) = \y a Ja,M a ) have 
the same parity. In this case, the g — > a transition is 
induced by the dc electric field [24] and the mixing of the 
states \g) and \a) with each other is characterized by the 
induced dipole matrix element 



(aid. 



dir. 



(8) 



We assume d m > since any complex phase can be incor- 
porated into the states \a) and \g). Note that the value 
of d- m depends on the dc field Ed c and the light polariza- 
tion e. In particular, d; n — > as Ed c — > 0. Therefore, the 
dynamic field does not cause mixing of the states \g) and 
| a) in the absence of the electric field. The polarizabil- 
ity a-g in Eq. ^ represents the ac polarizability of the 
unperturbed ground state \g), provided the effects of the 
dc field on the ac polarizability can be neglected. In this 
case, a-g = a g . 

In general, ac polarizabilities depend on the polariza- 
tion e and frequency u> of the external light field. The 
polarizability of an arbitrary atomic state \£) can be de- 
composed into three terms: 

M , 
ag = ao (7 J) + « ai (7 J) — (ex e ) • z 

. 3M 2 - J(J+1) 3|e-z| 2 -l 
+ " 2(7J) J(2J-1) 2 ' (9) 

The quantities aoiyJ), 0:1(7 J), and 0:2(7"/) are referred 
to as the scalar, vector, and tensor polarizabilities, re- 
spectively [3]. The scalar, vector, and tensor polarizabil- 
ities are independent of the magnetic quantum number 
M and the polarization e, and hence are independent of 
the choice of quantization axis and field geometry. How- 
ever, they depend on the light frequency w, as described 
in Appendix |A| 

Hereafter, we abandon the use of the overline to distin- 
guish between perturbed and unperturbed atomic states. 
Despite the lack of an overline, quantum states should be 



FIG. 2: Parameters of atomic trajectory. The atom's position 
r(i) and velocity v are shown for some time t > and have 
been projected onto two planes: (a) the plane formed by the 
vectors k and v, where k is the wave-vector of the standing 
wave of light, and (b) the plane normal to k and containing 
the origin O. A contour plot of the spatial extent of the light 
intensity is shown for reference; white and purple indicate 
minimal (zero) and maximal intensity, respectively. Here roy 
and ro± correspond to the components of r(t) that are parallel 
and perpendicular to k when t — 0. 



interpreted as atomic states that have been perturbed by 
the static electric and magnetic fields Ed c and B, unless 
otherwise noted. 

We consider a system in which atoms in the excited 
state undergo spontaneous decay to the lower states \b) 
and |c) with corresponding rates r a b and T ac — T a — r a t, 
where T a is the natural linewidth of the state \a). A 
schematic of the relevant energy level structure is shown 
in Fig. [l] As atoms decay from \b) down to \g), they 
emit fluorescent light of frequency uj^g. The LSM in- 
volves both the simulation and measurement of the spec- 
tral lineshape of the g — > a transition. In this context, 
the "spectral lineshape" refers to the probability of emis- 
sion of fluorescent light of frequency uibg as a function 
of laser frequency w. Although polarizabilities a g and 
a a depend on to, we assume that they are effectively con- 



stant for 



The LSM can be applied to any atomic 



system with the energy level structure shown in Fig. [T] 



III. SPECTRAL LINESHAPE 

The ac Stark shifts cause the resonant frequency of 
the g — > a transition to shift as atoms travel through 
the standing wave. As a result, the spectral lineshape 
depends heavily on the details of the light field. Because 
the light field amplitude E(r) is not spatially uniform 
[see Eq. ([2| and the discussion thereafter] , an atom with 
coordinate r(t) will experience a time-dependent electric 
field in its rest frame. 

Assuming constant velocity, the atom's position is 
r(t) = vt + Tq where v is the atomic velocity and Tq 
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is the position of the atom when t = 0. The origin O is 
chosen to be at one of the nodes of the standing wave. A 
diagram of the geometry is shown in Fig. [2] The time- 
dependent field experienced by the atom is given by 

E[r(t)} = E cos(kvnt + kr n) e-^ t2+r ^ )/r ^ (10) 



where E(r) is given by Eq. ([2|. Here v\\ — k • v and 

v± = | v— k vu I are the components of the velocity that are 
perpendicular and parallel to k. Similarly, r || and ro_i_ 
correspond to the perpendicular and parallel components 
of the position ro. 

The total time dependence of the electric field in 
Eq. ([I]) is due to the fast oscillation of the light at fre- 
quency uj and the slow modulation of the amplitude with 
a frequency kv\\ . The amplitude modulation is addition- 
ally characterized by a temporal Gaussian envelope with 
a characteristic width /v± , which is the amount of time 
an atom spends within the radius of the standing wave. 
We consider non-relativistic atoms for which the condi- 
tions uj 3> kv\\ and u v±/rh are valid, and the optical 
oscillations are much faster than the modulation of the 
amplitude E[r(t)]. In this case, the ac Stark shift is ob- 
tained by substituting Eq. (10 1 into Eq. 0. 



A. Absorption profile 

In order to gain a qualitative understanding of the 
physics, we make the simplifying assumption 



E(t) fa £; cos(fc7;||t). 



(11) 



In this case, the Stark shifts of the ground and excited 
states lead to the following shift of the resonant energy 
of the g — > a transition: 



£ ac (a) ~ £ac(flO 



1 



where 



i ag E 2 cos 2 (kv ll t), (12) 



(13) 



is the difference of the polarizabilities of the ground and 
excited states. From an atom's perspective, this is equiv- 
alent to a polarizability-dependent frequency modulation 
of the two counterpropagting light fields. Thus, the fea- 
tures of the g — > a lineshape can be understood by study- 
ing a related system: stationary atoms with fixed en- 
ergy levels in the presence of two counter-propagating, 
frequency-modulated electric fields. In this subsection, 
we turn our attention to such a system. 

The frequency-modulated electric fields have instanta- 
neous frequencies ui + and ui- given by 



uj± (t) — uj ± fcu|| + u! m A m cos(uj m t), (14) 



where 



a a E 2 

uj m = 2kvu and A m = ag , 
" onkvw 



(15) 



are the modulation frequency and modulation index, 
respectively. Equation (14) includes the term ±kv\\ 
which accounts for the Doppler shifts of the frequen- 
cies of the two counter-propagating waves. To derive 
Eq. ( |14[ ), we used a trigonometric identity to write 
cos 2 (kv\\t) = (1/2)[1 + cos(2kv\\t)] and we neglected the 
time-independent term because it can be interpreted as 
an overall shift of the optical frequency: u — > uj + uj m A m . 
The instantaneous frequency ui± (t) is characteristic of a 
light field with a time-dependent phase [25] . Such a field 
is given by 



E, 



(l/2)E e i ^ u}±kvi1 ) t+j4m sin ("n,t)] ^ 



(16) 



where the factor of 1/2 is included so that the total field, 
which is the sum of two traveling waves, has an amplitude 
of Eq. The effective field E± can be decomposed in the 
following way: 



E±= J2 (V2)£b-/»(4n)e 



i(o;±fcu|| +nu; m )i 



(17) 



where J n are Bessel functions of the first kind. Thus 
the effective field consists of a principal field (ri = 0) 
which oscillates at a frequency ui ± kv» , and infinitely 
many sidebands (n ^ 0) which oscillate at frequencies 
uj ± + nuj nl . The amplitude of the electric field of 
nth sideband is (l/2)E J n (A m ). 

The total electric field E to t seen by the atom is the 
sum of the two counter-propagating light fields. To add 
the fields, the summation index in Eq. (17) is changed 
from n to n + 1 and the field _E_ is expressed as 



E_ = 



OO 

E 



(l/2)E J n+1 (A 



where —kv\\ + uj m = +kv\\ because the modulation fre- 
quency uj m is exactly twice the Doppler shift. Hence the 
sidebands of the two counter-propagating waves overlap 
and the total field is given by 



E u 



(18) 



71 — — OO 



where 



E 

E n = — [J n +i + J n ] and ui n = ui + kv\\ +2n kv\\ . (19) 

In particular, the first-order sidebands from one field cor- 
respond with the carrier of the other [23] , resulting in an 
absorption profile with a polarizability-dependent distor- 
tion. The "absorption profile" is a plot of the transition 
rate as a function of uj. 

The rate of the g —¥ a transition is given by 



K 2 



T a /2ir 

(wn - u ag y + (r a /2y 



(20) 
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than the peak on the right. The following conditions need 
to be met in order for the distortion to be observed: 



-3-1 1 3 5 

u - uj ag [uj-d] 

FIG. 3: Absorption profile. Shown is a plot of the transition 
rate 1Z as a function of the laser frequency u) for the case of a 
single atom with parallel speed vu = u>u/k, where ojo is the 
Doppler broadening of the line. The transition rate is normal- 
ized by IZo to have a maximum value of unity. Here the mod- 
ulation index satisfies Am < 0. For Am > 0, the peak on the 
left is taller than the peak on the right. In the insert, the av- 
erage transition rate 1Z is shown in three cases: (a) the condi- 
tion r a /4 <C ijJd is violated, (b) the conditi on Fi uj-p < |a as |_Eo 
is violated, and (c) both conditions in Eq. (22 1 are satisfied. 



To generate these curves, we made use of the approximation 
k ~ k ag which is valid in the near-resonant regime (uj ~ u ag ). 



Equation (20) is valid in the weak excitation limit, that 
is, when the excitation rate 1Z is much smaller than all 
other relevant rates. To derive Eq. (20), we neglected 
the interference of different harmonic components, e.g., 
£ n e 4 ""* and 8 n ie %Unlt . Such terms contribute small cor- 
rections to the transition rate which do not affect the 
qualitative behavior of the absorption profile. A plot of 
the absorption profile is given in Fig. [3] The single-atom 
absorption profile is clearly asymmetric about the atomic 
resonance (uj — uj ag )- 

For an ensemble of atoms, the absorption profile is ob- 
tained by averaging the transition rate ( 20 ) over the ve- 
locity distribution. The average rate is 



K : 



Tl(v\\) wi,(v\\,vc)dv\\ 



(21) 



where wi,(vn,vc) is the appropriate probability distri- 
bution for the parallel velocity v\\ , and vq is a charac- 
teristic speed in the k direction. A plot of the average 
transition rate is shown in the insert of Fig. [3] For the 
insert, w^lvu , vc) is taken to be a Lorentzian distribution 
with a full-width at half the maximum value (FWHM) of 
vc = uJn/k, where wd is the overall Doppler broadening 
of the line. The resulting absorption profile is similar to 
experimentally observed lineshapes in Yb (Section IV). 

The absorption profile exhibits a polarizability- 
dependent feature: a dip that separates the profile into 
two distinct peaks. The sign of the polarizability a ag de- 
termines whether the peak on the left is larger or smaller 



Ta/4 <C wd and hw^ < 



a ag \E . 



(22) 



The former condition ensures that the sidebands of the 
FM waves can be resolved. The latter ensures that the 
amplitude of the first-order sidebands is not negligible 
compared to the amplitude of the carrier, that is, A m =/= 
for W|| ss vq- If cither wd S Ta/4 or hu>r> ^> |a ag |-Eg, then 
the asymmetric distortion will be suppressed, as can be 
seen in the insert of Fig. [3| In this case, the methods de- 
scribed here cannot be used to measure the polarizability 
a ag . However, the LSM can still be used to measure a ag 
in the absence of the distortion by comparing displace- 
ments of the central peak of the lineshape. 

By omitting the Gaussian envelope exp(— v\t 2 jrl^) in 



Eq. (11), we neglect effects of the atom's finite tran- 
sit time, such as broadening of the spectral line |26j . 
Nonetheless, the transit time must satisfy the following 
restrictions: 



v±/r h < d hl E Q /h and uj_/r b <WD, 



(23) 



where di n Eo/h is the Rabi frequency of the g — > a transi- 
tion. The former condition represents a system in which 
atoms have enough time to undergo excitation to the up- 
per state \a), as will be discussed in Section MB The 
latter condition ensures that the Doppler broadening is 
sufficiently large that most atoms travel through many 
nodes and antinodes of the standing wave during their 
transit. 

Although the absorption profile provides a satisfactory 
illustration of the physics, it cannot be used to measure 
the ac polarizabilities. The transition rate presented in 
Eq. ( 20 1 is valid only in the weak excitation limit and 
therefore cannot account for saturation effects. More- 
over, Eq. ( 20 ) does not properly take into account in- 
terference of different probability amplitudes, nor does it 
include finite transit-time effects. A more complete pic- 
ture is required to generate a theoretical lineshape that 
can be fitted to experimental data. Such a picture is 
achieved by the following model. 



B. Fluorescence probability 

Hereafter, we return our attention to the original sys- 
tem: moving atoms illuminated by light with a fixed fre- 
quency uj. The spectral lineshape of the g — ¥ a transi- 
tion is modeled by computing the probability of emis- 
sion of fluorescent light of frequency uib g as a function 
of laser frequency ui. The computation involves three 
steps. First, the time-dependent population of the 
state | b) is computed by numerically solving the optical 
Bloch equations (OBE) for the case of atoms traveling 
through a non- uniform field (see Fig. [2]). Second, the 
probability of fluorescence J- is determined by integrat- 
ing the time-dependent decay rate T^ptb with respect to 
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time. Finally, the average fluorescence probability T is 
computed by taking a weighted average of JF with respect 
to the atomic velocity v and the offset r . 

Although the LSM can be used in conjunction with 
any atomic source, our model makes use of distributions 
that are appropriate for a collimated beam of thermal 
atoms traveling in a direction that is orthogonal to the 
standing wave. In this case, v± represents the compo- 
nent of the atom's velocity along the atomic beam, and 
V\\/v± represents the angular spread of the beam. The 
corresponding velocity distribution is 



In terms of the dimcnsionless parameters, the condi- 



tions presented in expressions (22) reduce to S 2 3> 1/2 
and I Si I > 1/8. When either of these conditions is vi- 
olated, the characteristic dip in the lineshape is sup- 
pressed, as can be seen in Fig. [4j Likewise, condi- 
tions p3]) reduce to S 3 < y/So~ and S 3 < (l/ 2 )^- 



w(v) = ■wt(v±,v t )wi j (v\i,vc), 



where 



w t {v_l,v t ) = 2(v ± /v T )e 



(24) 



(25) 



is the distribution of velocities appropriate for thermal 
atoms escaping from a hole, and mi^Wuiic) is the ve- 
locity distribution appropriate for a collimated atomic 
beam. Here i>t = ^/^ksT/M is the thermal speed of the 
atom, T is the temperature of the oven, M is the atom's 
mass, Icb is Boltzmann's constant, and vq is the charac- 
teristic speed determined by the atomic-beam collimator. 
To model the effects of a vane collimator, we approximate 
the spread of parallel velocities by a Lorentzian distribu- 
tion with a FWHM of vc = /k, where ojd is the overall 
Dopplcr broadening of the line. 

In the following model, we use dimensionless param- 
eters. Dimensionless parameters ease computation, and 
potentially facilitate the application of the model to sev- 
eral different atomic systems. Throughout, we make 



Whereas the absorption profile discussed in Section [ill A| 
was valid only in the weak excitation limit (So ^ 1) 5 the 
model of the fluorescence can accommodate large satu- 
ration parameters. 

Let p nm be the elements of the density matrix in the 
atom's rest frame for states n, m £ {g, a, 6, c}. We as- 
sume that the rotating wave approximation holds and 
dynamic interactions between states other than \g) and 
I a) can be neglected. In this case, the dimensionless op- 
tical Bloch equations (OBE) for the configuration shown 
in Fig. [l]are [27] 



Paa 

Pag 
Pbb 



--fl(Pag - Pga) - Paa, 



-^(pgg - p aa ) - 

-GbPbb + G ao p aa , 
~G c p cc -\- G a cPaai 



(l-2iA) Pag , 



(30a) 

(30b) 

(30c) 
(30d) 



where p nm = d{p nm ) / dr . The remaining density matrix 



elements p gg and p ga 
and p ga = pig. Here 



are determined from J2 n Pnn = 1 



n = 



1 



\d-e\g)E(t) 



= V^b/(C,u), 



(31) 



k aq which is valid in the near- is the Rabi frequency, 



J "9 



the approximation k 
resonant regime (u) ~ 

Time is measured in units l/r a . We define the di- 
mensionless time t = T a t and decay rates Gb = Tb/T a , 
G c = r c /r a , Gab = Tab/Ta, and G ac = 1 - G ab - We fur- 
ther define the dimensionless perpendicular and parallel 
velocities u± = v±/v^ and u\\ = v\\/vc, and the dimen- 
sionless perpendicular and parallel offsets ^ = r j_/rb and 
ip = fcroy, respectively. 

we in- 



A = 



r„ 



8(a) -£{g) 



= 6 + 2S 1 S 2 f{£,,u)\ (32) 



is the detuning of the laser light from the resonance, 
8 = {lu — Lu a g)/T ai and the function /(£, u) is defined by 



/(€,u) = cos[($ 2 /2)tt||T + de 



-HS 3 u ± r) 2 +e} 



(33) 



Consistent with the discussion in Section III A 



troduce the following dimcnsionless parameters: the sat- 
uration parameter So, characteristic modulation index 
5*1, and Doppler parameter S2, defined by 



where £ = (<p,£) and u = (ui\,u±_). We further assume 
that all atoms initially occupy the ground state. 

The probability that an atom will emit fluorescent light 
of frequency LObg in a time interval [r;,Tf] is given by 



5 = [d in E /(hT a )} , 

51 = a ag E^/(8fuJr>), 



(26) 
(27) 



G h pbb{r) dr, 



(34) 



and 



S 2 = 2w D /r a , (28) 

respectively. We define an additional parameter S3 by 

S 3 = v T /(r h T a ). (29) 

Note that r^,/v± is the time that an atom spends within 
the radius of the light field and hence l/(S 3 u±) repre- 
sents the dimensionless transit time. 



where Gb Pbb{T~) is the time-dependent rate of the b — > g 
decay. The choice of integration interval [71, Tf] depends 
on both the transit time and the characteristic time of the 
fluorescent decay after the atoms leave the light field. In 
the case where the a — > b — > g decay time is shorter than 
the transit time, it is appropriate to define the integration 
limits by — t\ = T{ = 3/(S 3 u±). The factor of 3 ensures 
that the atom is "far" from the standing wave at the 
integration limits. To model a system with slower decays, 
the integration limit Tf must be extended. 
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FIG. 4: Results of simulation. The simulated lineshape J- (5) 
is shown for two different values each of the modulation index 
Si and the Doppler parameter S2, as indicated by the insert: 
(a) the condition S2 ^> 1/2 is violated, (b) the condition 
Si I > 1/8 is violated, and (c) both conditions are satisfied. 
For each curve, the saturation parameter is So = 1.0. The 
lineshapes are normalized by the maximum peak height Fq. 
To ease comparison, the middle and top curves are shifted 
vertically by 0.7 and 1.0, respectively. 



The fluorescence probability depends on the parame- 
ters of the atomic trajectory, that is, J- = J-(£,u). We 
define the average probability of fluorescence by 



T 



J-(£, u) w(£, u) e?£ du, 



(35) 



where w(£, u) is the probability distribution associated 
with the atom's initial position and velocity. In our 
model, we assume w(£, u) = itfi(£) 1^2(11), where the 
distribution w\ is a uniform distribution over the inter- 
vals £ G [—3,3] and ip e [— tt,it]. The finite integra- 
tion limits are justified by the following properties of 
the system: First, the amplitude of the standing wave 
drops to less than 0.01% of its maximum value when 
|£| > 3. Therefore, atoms will only pass through the 
light if |£| < 3. Second, ip constitutes a phase shift of the 
electric field which is unique only for ip G [— 7r,7r]. Con- 
sistent with Eqs. (24 1 and (25), the velocity distribution 
satisfies u^u) = Wt(u±, l)tujii(ui|, 1), where wt(u±,1) is 
the velocity distribution for atoms escaping from a hole 
with unit thermal speed, and wi,(u\\,l) is a Lorentzian 
distribution with FWHM of 1. 

The fluorescence probability is a function of the detun- 
ing 5 from the atomic resonance, that is, 



(36) 



where S = (So, Si, S2, S3) and G = (G a b, Gb, G c ) are pa- 
rameters. We refer to a plot of J 7 (5; S, G) as a function 
of 6 as the "simulated lineshape" of the g — > a transition. 
Three such plots are shown in Fig. [4] In Fig. |4j curves (a) 
and (b) demonstrate the suppression of the polarizability- 
dependent distortion when either the Doppler parameter 



52 or the modulation index 5i is too small. The distor- 
tion is most pronounced in curve (c), for which 52 > 1/2 
and 15*11 > 1/8. The simulated lineshapes in Fig. [3] are 
qualitatively similar to the absorption profiles in the in- 
sert of Fig. [3j as expected. 



C. Numerical procedure 

The numerical procedures described here are valid for a 
variety of atomic species. However, the simulations were 
performed with parameter values appropriate for the Yb 
system described in Section [iVj 

We used a stiffly stable Rosenback method [35J to nu- 
merically solve a system of equations related to Eqs. (30) 

3/(53 u ±)- This system of 
The Rosenback 



and (34), with — t; = n 
equations is described in Appendix [B] 
method involves two tolerances-denoted atol and rtol 
in Ref. 1281 -which were both set to 10~ 5 . The multi- 



dimensional integral in Eq. ( 35 1 was computed using an 
adaptive Monte Carlo routine |28j. In our implementa- 
tion, the integration routine involves 10 5 evaluations of 
the integrand. For various values of S, the average esti- 
mated error was less than 1% of the value of the integral. 

For computational purposes, we restricted the inte- 
gration to the following finite domain: £ G [—3,3], 
p € [— 7r,7r], u± G [0,3], and uii <E [—6,6]. The subdo- 
mains for £ and (p were discussed after Eq. ( 35 1 . The fi- 
nite integration subdomains for un and u± are justified as 
follows: Atoms with a large parallel speed \u» | experience 
a Doppler shift that is much larger than the characteris- 
tic Doppler broadening of the spectral line. Such atoms 
only contribute to the wings of the lineshape, where \S\ 
is large and the probability of fluorescence is very small. 
Moreover, for a Lorentzian velocity distribution with unit 
FWHM, |u||| < 6 for about 95% of atoms. 

On the other hand, atoms with a perpendicular speed 
that satisfies u± 3> v / 5o/5*3 are moving so fast that the 
transit time is much smaller than the inverse Rabi fre- 
quency. Such atoms do not spend enough time in the 
light field for the g — > a transition to be realized. Since 
most atoms travel at or near the thermal speed u± = 1, 
the condition ^/So/S^ > 1 represents a system in which 
most atoms have enough time to interact with the light. 
We assume that atoms with speed u± > 3 do not con- 
tribute significantly to the lineshape. Note that u± > 3 
for only about 0.1% of atoms. Finally, we ignore counter- 
flow of atoms in the atomic beam by requiring u± > 0. 

For fixed G = (0.35,0.45,0) and 5 3 = 0.53, we com- 
puted the average fluorescence J-{8; S, G) for 100 dis- 
crete values of 5 G [5_,5+], and various discrete values 
of S G [0.1,10], 5i G [0.01,1], and 5 2 G [1,100]. Here 
6± = 25i5 2 ± 3[(l/2)5* 2 +G b + 1}. The results were in- 
terpolated using cubic splines to approximate the contin- 
uous function T(S;S,G). Three curves which are typi- 
cal of those produced by this procedure are presented in 
Fig. [4] 

The LSM involves fitting the simulated curve 



J- (5; S, G) to the observed lineshape to determine best- 
fit values of So, Si, and S2. From the best-fit values, the 
following three quantities can be calculated: the polar- 
izability difference a ag , the circulating power P of the 
standing wave 1 , and the Doppler broadening wd of the 
g — > a transition which are given by 



P = 



1 

Ic 



x SiS 2 /S 



x S , 



and 



lju = (r o /2) x s 2 , 



(37) 
(38) 

(39) 



respectively. Here eo is the permittivity of free space. 
The present implementation of the LSM uses Mathemat- 
ical nonlinear regression routine to determine the best- 
fit values of S , Si, and S 2 . Alternatively, Eq. (38 1 can 
be solved for d[ n in terms of P and So- Thus the LSM 
can also be used to measure the induced dipole moment 
when the power is known, as was done in the previous 
application of the LSM [25] , 

The treatment of systematic errors and statistical un- 
certainties is straightforward. In Eq. (37 1, for instance, 
the uncertainty in the dimensionful quantity df n /(HT a ) is 
due solely to systematic effects, whereas the uncertainty 
in the term S1S2/S0 arises from statistical uncertainties 
in both the observed signal and the fitting algorithm. 
The total uncertainty of the quantity a ag is obtained by 
adding these independent uncertainties in quadrature. If 
the signal-to-noise ratio (SNR) of the observed lineshape 
is sufficiently high, then the error of the measurement of 
a ag will be dominated by the uncertainties of the known 
quantities d- ln and T a . 



IV. APPLICATION TO YTTERBIUM 

The electronic structure of Yb is shown in Fig. [5j The 
low-lying energy eigenstates of Yb match the structure 
shown in Fig. nlunder the following mapping: \g) = |'So), 
I a) = | 3 Di), 1 5) = | 3 Pi), and |c) = | 3 P ). Therefore, the 
LSM can be used to measure the difference in ac po- 
larizabilities of the upper state | 3 Di) and the ground 
state |'S ) at 408 nm by analyzing the lineshape of the 



1 The standing wave is formed by two counter-propagating waves 
of light. The circulating power of the standing wave is de- 
fined as the average power of a single traveling wave. The 
electric field of the wave propagating in the ±k direction 
is given by Ei(r, t) = (1/2) Eo exp[— (r±/r h ) 2 ] cos(fcr|| ± uit). 
The corresponding time-averaged intensity in SI units is 
I(r±) = (l/8)ceoE 2 exp[— 2(rj_/r^) 2 }, where the time aver- 
age is taken over a single period of oscillation. There- 
fore, the circulating power of the standing wave is given by 
P = / °° I(r±) 2irr ± dr ± = (1/16)^ ce E 2 . 
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FIG. 5: Ytterbium energy levels. Shown are the low-lying en- 
ergy eigenstates of Yb and the electronic transitions relevant 
to the polarizability experiment. 



408-nm 'So 3 Di transition. In this case, the line- 
shape is measured by observing the 556-nm fluorescence 
of the 3 Pi — > 1 Sq decay. The intermediate state | 3 Po) is 
metastable and hence T c « 0. 

There is an additional decay of | 3 Di) to the metastable 
state | 3 P2), which is not shown in Fig. [I] The state |c) 
can represent multiple metastable states, including both 
| 3 P i ) and | 3 P2) . In this interpretation, r ac is the rate of 
decay of | 3 Di) to all metastable states. 

The highly forbidden 'So — > 3 Di transition is in- 
duced by the Stark mixing technique. This technique 
involves the application of a static, uniform electric field 
Edc which mixes the upper state | 3 Di) with opposite- 
parity states, predominantly the |'Pi) state. The effec- 
tive dipole moment d m associated with the Stark-induced 
transition is given by |23j 



Pag | (Edc 



>-M„ 



(40) 



where j3 ag = 2.18(10) x 10~ 8 e • a • (V/cm) -1 is the vec- 
tor transition polarizability of the 'So — > 3 Di transition, 
e is the elementary charge, et is the Bohr radius, and 
(Edc x e)W is a spherical tensor of rank one 2 . Here M a 
is the magnetic quantum number of the | 3 Di) state. 

Because the angular momentum of the ground state 
| 'So) is J g = 0, only the scalar term in Eq. ^ 
contributes to the polarizability of I 1 So)- That is, 
a g = Q ; o( 1 So) and so a g is independent of the geome- 
try of the applied fields. Hence the dependence of a ag on 
the field geometry is due entirely to the vector and tensor 
polarizabilities ai( 3 Di) and «2( 3 Di) of the excited state 
| 3 Di). In this case, the LSM is sensitive to the difference 
of the scalar polarizabilities of the ground and excited 



2 Let AW be the rank-one spherical tensor associated with the 
Cartesian vector A. Then the components of AW are given by 



d 



A, and A 



(i) 
±1 



± iA y )/y/2. 
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states. However, the vector and tensor polarizabilities 
of the excited state can be measured unambiguously by 
varying the polarization of the standing wave. 

A recent calculation of the polarizability of the 
ground state I 1 So) at 408 nm yielded 3 

aoCSo) = 0.218(11) Hz • (V/crnT 2 . (41) 

Calculations of the polarizability | 3 Di) at 408 nm are 
complicated by the potential existence of odd-parity 
eigenstates with energy close to twice the energy of a 
408-nm photon. Such states could lead to a resonantly 
enhanced polarizability of the 3 Di state. The energy 
spectrum in this region (which is below the ionization 
limit) is very dense due to the excitation of 4/ orbitals. 
The knowledge of the energy spectrum is far from com- 
plete in this region. This provides one of the motivations 
for determining the polarizabilities experimentally. 

The first implementation of the LSM [23] was used to 
measure the quantity 

a\ g = -0.312(34) Hz ■ (V/cm)- 2 , (42) 

were a l ag — a2( 3 Di) + a ( 3 Di) — «o( 1 So). Here the su- 
perscript "I" is introduced to distinguish this measure- 
ment from the results of the present work. To determine 
the tensor contribution «2( 3 Di) unambiguously requires 
a second measurement of a different combination of scalar 
and tensor polarizabilities. This is accomplished in the 
present work by the application of a dc electric field that 
is parallel to the standing wave (Fig. [6| , whereas the pre- 
vious measurement was performed with a dc field that 
was perpendicular to the standing wave. In addition, the 
current experiment includes a strong magnetic field not 
present in the previous case. The magnetic field makes 
possible the measurement of the vector polarizability, as 
discussed in Section IVl 

A. Experimental apparatus and field geometry 

The details of the experimental apparatus were re- 
ported elsewhere [14] , and only a brief description is pro- 
vided here. A schematic of the setup is shown in Fig. [6] 
A beam of Yb atoms is produced by a stainless-steel oven 
loaded with Yb metal, operating at 500 °C. The oven is 
outfitted with a multislit nozzle, and there is an exter- 
nal vane collimator reducing the spread of the atomic 
beam in the z-direction. Downstream from the colli- 
mator, atoms enter a region with three external fields: 
a uniform, static magnetic field B; a uniform, static 
electric field E^; and a non- uniform, dynamic electric 
field E(r,i). 



Equation (|7jl implies that ct ag has units of energy per squared 
electric field. However, in this work a ag is normalized by 2irh 
and presented in units of frequency per squared electric field. For 
a more thorough discussion of unit conventions, see Ref. [5]. 



y 




FIG. 6: Experimental apparatus. A collimated beam of ytter- 
bium atoms interacts with a standing wave of light in the pres- 
ence of dc electric and magnetic fields. The light is resonant 
with the 408-nm 1 So — > 3 Di transition, and the 556-nm fluo- 
rescence is detected by the photomultiplier tube (PMT). The 
neodymium (Nd) magnets are axially magnetized in the z- 
direction. The standing wave is generated in a power buildup 
cavity (PBC). With the exception of the PMT, the apparatus 
is contained in a vacuum chamber. 



The magnetic field B is generated by a pair of axially 
magnetized neodymium (Nd) magnets. These magnets 
produce a field with sufficient strength (more than 50 G) 
to completely isolate the Zeeman sublevels of the upper 
state | 3 Di). The electric field Ed c is generated by two 
wire-frame electrodes separated by 2 cm. The ac electric 
field E(r, t) is due to standing-wave light at the transition 
wavelength of 408.346 nm in vacuum, which is produced 
by doubling the frequency of the output of a Tksapphire 
laser (Coherent 899). About 7 mW of 408-nm light is 
coupled into a power buildup cavity (PBC) with finesse of 
approximately 15,000. The PBC is an asymmetric cavity 
with a flat input mirror and a curved back mirror with a 
50-cm radius of curvature. The separation between the 
mirrors is 22 cm. 

Fluorescent light with a wavelength of 556 nm is col- 
lected with a light guide and detected with a photomul- 
tiplier tube (PMT). With the exception of the PMT, the 
entire apparatus is contained within a vacuum chamber 
with a residual gas pressure of 3 x 10~ 6 Torr. 

As can be seen in Fig. [6] the fields B and Ed c point in 
the z direction. Likewise, the standing wave is oriented 
along the z-axis. The light field E(r,t) lies in the xy- 
plane. For this geometry, the transition to the upper 
state | 3 Di;M a = 0) is suppressed. 

The polarization of the light field is of the form 
e = e x x + e y y, where \e x \ 2 + \ e y\ 2 = 1- To further char- 
acterize the polarization, we introduce three parameters: 
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FIG. 7: Results of fitting program. Shown is a comparison of 
the observed and the simulated lineshapes of the 1 So — > 3 Di 
transition. The data correspond to Run 3 in Table [I] Also 
shown are the residuals, magnified by a factor of 10. The 
residuals are the difference of the data and the simulation. 



polarization angle 9, degree of ellipticity 4>, and handed- 
ness h, which are given by |30j 



and 



tzn20 = 2Re(e x e* y )/{\e x \ 2 ~ \e y \ 2 ), 
sin20 = 2|Im(e ;E e;)|Hi(exe*)-z| 



-sgn[t(e x e*) ■ z]. 



(43) 
(44) 

(45) 



Linearly, circularly, and elliptically polarized light are de- 
scribed by 4> = 0, 4> = 7r/4, and < (p < 7r/4, respec- 
tively. The sense of rotation is determined by h: left 
and right-handed polarizations correspond to h = +1 
and h = —1, respectively. Substituting Ed c = E Ac z and 
e = e x x + t y y into Eqs. ( 40 ) and Q yields 



din = PagE Ac y/{\/2) [1 + hM a S hl(20)] , (46) 



and 



a a = a ( 3 Di) - / l M a sin(20)ai( 3 D 1 ) - -a 2 ( 3 Di), (47) 



for M a = ±1. Here we have used Eqs. (44) and (45) to 
eliminate e x and e y in favor of the degree of ellipticity (f> 
and the handedness h. For this geometry, both d- m and 
a a are independent of the polarization angle 6. 

According to the geometry in Fig. [6j the component 
of the atom's velocity that is perpendicular to the stand- 
ing wave is v± = (v 2 + v 2 ) 1 / 2 . The output of the oven 
is about 6 mm in the y-direction, and is located more 
than 20 cm away from the standing wave. In order for 
an atom to pass through the standing wave, its veloc- 
ity components must satisfy \v y /v x \ < 0.015. Therefore, 
the approximation v± — v x is valid and the use of the 
thermal distribution given in Eq. (25) is justified. How- 



ever, for a high-precision measurement, the effect of this 
approximation needs to be investigated. 



TABLE I: Summary of results. The uncertainty in the fitted 
parameters is the asymptotic standard error returned by the 
fitting algorithm. 



Run 


co-10 2 


ci-10 2 


So 


Si 


s 2 


S1S2/S0 


1 

2 
3 
4 


5.1(1) 
5.0(1) 
5.4(1) 
5.0(1) 


5.5(4) 
3.2(3) 
7.0(3) 
4.7(3) 


2.83(17) 
2.76(16) 
2.68(12) 
2.56(17) 


-0.208(3) 
-0.215(3) 
-0.187(2) 
-0.195(3) 


34(1.1) 
30(1.1) 
34(1.0) 
29(1.1) 


-2.53(16) 
-2.33(17) 
-2.35(13) 
-2.22(17) 


Avg.: 






2.70(7) 


-0.199(1) 


31.8(5) 


-2.36(8) 



B. Data analysis 

The data were acquired over four separate experiments 
(runs). During each run, approximately 2000 lineshapes 
were recorded at a rate of about 100 ms per lineshape 
and an average lineshape was computed. This procedure 
resulted in 4 lineshapes, each with an effective integration 
time of about 200 s. A typical lineshape is shown in 
Fig. [7] After fitting the theoretical model to the data, the 
statistical uncertainties of each run were scaled to give a 
reduced \ 2 01 unity. The scale factor varied between 4.0 
and 4.7 depending on the run. The resulting error bars 
are shown in the figure. The relatively large scale factor 
indicates that the accuracy of the fit is dominated by 
either systematic distortion of the profile during the scan, 
or profile features that are neglected in the theoretical 
model, but not by the statistical uncertainty of the signal. 

In the present experiment, the oven temperature was 
T = 500(50) °C and the magnitude of the dc electric 
field was E^c — 4.24(6) kV/cm. The atoms intersect the 
standing wave in the middle of the PBC where the radius 
of the light beam is r^ — 196(5) /im. We used linearly 
polarized light with cp = 0(1)°. Then Eqs. ([37]) through 
(39) become 



= [0.0668(86) Hz • (V/cm)" 
P = [5.02(98) W] x S , 



x S 1 S 3 /S Q , (48) 
(49) 



and 



cj d = [2tt x 0.209(17) MHz] x S 2 , 



(50) 
Here 



where c£ g = a ( 3 Vi) - (l/2)a 2 ( 3 D 1 ) - c^So) 
the superscript "II" is used to distinguish the re- 
sults of the present work from the previous measure- 
ment. The dimensionless decay rates are G a b — 0.35(4), 
Gi, = 0.45(4), and G c = 0, and the parameter S3 is given 
by S 3 = 0.53(5). 

As part of our analysis, we normalized the observed 
and simulated lineshapes by their maximum values. In 
addition, we introduced two calibration parameters to 
the simulated curve: 

T{x\ S, G, c) = (1 + c ) JX(1 - c{)x; S, G) - c , (51) 

where c = (co,ci). Here cq accounts for the background 
of the observed signal. The parameter ci is a scaling fac- 
tor that accounts for any variation in the calibration of 
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TABLE II: Error budget. Shown are the factors contributing 
to the uncertainty of the measured value of a 1 ^ . 



Factor 


Uncertainty (%) 


Vector transition polarizability (j3 ag ) 


9 


Lifetime of 3 Di 


8 


Light polarization 


4 


DC electric field 


3 


Simulation and data fit 


3 


Total (in quadrature) 


13 



the frequency axis of the data relative to the simulation. 
Such deviations could arise due to misalignment of the 
atomic beam relative to the axis of the PBC, deviations 
in the perpendicular velocity distribution, or the uncer- 
tainty of the timescale 1/T a . In practice, both c$ and c\ 
represent very small corrections, with typical values on 
the order of 0.05 (see Table [1. 

The results of the fitting for each run are given in Ta- 
ble |TJ Combing these results with Eqs. (48) through (50) 
yields 



= -0.158(21) Hz • (V/cm)" 



(52) 



P = 13.6(2.6) W, and cj d = 2tt x 6.58(53) MHz. The 
measured values of the circulating power P and the 
Doppler broadening wd are consistent with (and more 
precise than) direct measurements of these quantities. A 
comparison of the data to the fit is shown in Fig. [7J The 
quality of both the data and the fit are sufficiently high 
that the uncertainty of the measured value of a l ^ g is pri- 
marily due to the uncertainties in the vector transition 
polarizability f3 ag and the linewidth r a . A summary of 
the factors that contribute to the uncertainty of the mea- 
surement are shown in Table HT1 

The quantity a^ g = a ( 3 Di) - (l/2)a 2 ( 3 D!) - ao( x S ) 
is a combination of the ac scalar and tensor polarizabili- 
ties of the states 3 D! and 1 So- The tensor polarizability 
of 3 Di is determined unambiguously by comparing the 
present measurement of a l ^ g with the previous measure- 



ment of a J ag , which was discussed after Eq. (|42|). We 
find 



a 2 ( 3 Di) = -0.103(26) Hz • (V/cm)" 



(53) 



and 



a ( 3 Di) - ao^So) = -0.209(18) Hz ■ (V/cm)" 2 . (54) 

Finally, the scalar polarizability of 3 Di is isolated by sub- 
stituting the calculated value cvo( 1 So), given by Eq. (41 ), 
into Eq. ( p54| ): 

a ( 3 Di) = 0.009(21) Hz • (V/cm)" 2 . (55) 

With the present accuracy, the scalar polarizability of 
3 Di is consistent with zero. In the presence of linearly 
polarized light, the polarizability of 3 Di is dominated by 
the tensor polarizability. However, for light with arbi- 
trary polarization, a 1 ( 3 D 1 ) also plays a role. The vector 
polarizability is the subject of ongoing experiments. 



V. SUMMARY AND OUTLOOK 

This work is part of a continuing investigation of po- 
larizabilities in Yb. The ac scalar and tensor polariz- 
abilities of the excited 3 D! state in Yb were measured 
independently for the first time. Ongoing experiments 
are focused on measuring the ac vector polarizability, for 
which there is currently no experimental or theoretical 
data. 

To measure the vector polarizability, the 1 So — > 3 Di 
transition must be excited using circularly polarized 
light, as can be seen in Eq. (47 1. Such a measurement 
requires control over the ellipticity of the light. In the 
present experimental setup, only two Zeeman sublevels 
(M a = ±1) are excited. The degree of ellipticity can 
be measured by comparing the relative strengths of the 
transitions to different sublevels. For purely circularly 
polarized light, only one sublevel is excited. This condi- 
tion is ideal for measurement of the vector polarizability. 

In this paper, we presented the next generation of 
the Lineshape Simulation Method (LSM) for measuring 
combinations of polarizabilities of the ground and ex- 
cited states in atoms. The LSM was originally developed 
specifically for Yb, but we have generalized the method 
for an arbitrary atomic system with the level structure 
shown in Fig. [7J For example, the LSM could be used to 
measure the polarizabilities of the 6S and 7S states in ce- 
sium by observing the lineshape of the 539-nm 6S — > 7S 
transition driven by a standing wave of light [11] . 
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Appendix A: Frequency dependence of dynamic 
polarizabilities 

The ac polarizability ae of the state \£) = YyJM) is 
given by Eq. ([9]). The scalar, vector, and tensor polariz- 
abilities depend on the light frequency w in the following 
way 011]: 



ai( 



1 


E 

t'^i 


\du>\ 2 


3h 


2J + 1 


l 


E 


\d u >\ 2 


2h 


2J + 1 


l 


E 


\d u >\ 2 


3h 


2J + 1 



Here the summation is over all states \£') — \j'J'M'} 
such that \l) and \C) have opposite parity. The functions 
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$Cb 3>i an d ^2 are given by 
$o(J, J') = <5jj< + <5j,j'+i + 



$i(J,J') = - 



J + 1 



-<y. 



and 



(2J-lR^ J(2J-l)^j, +1 

$2(J ' J) -^/Ti (j+i)(2j+3) 



(A4) 
(A5) 

(A6) 



where Sjji is the Kronecker delta. We emphasize that 
<I>i(J, J') differs from the expression found in Refs. [31 H] 
by a factor of J. The reason for this discrepancy is that 
we follow the convention for which the vector polariz- 
ability of a stretched state (M = J) is 01(7 J) instead of 
J x aii'jJ). 



( A3 ) are given by 



The frequency-dependent parts of Eqs. (|A1|) through 

(A7) 



1 



± 



1 



we>e — u) — iT(/ /2 uj£>£ + u> + iTfi /2 



Here dw = ('fJ\\d^\\'-f' J') is the reduced matrix el- 
ement of the electric dipole operator and d^ 1 ' is the 
spherical tensor associated with the electric dipole mo- 
ment. In the limit of large detuning from resonance 
(\w#t ±cj\ ^> T11), the functions Fpii reduce to 



F + - 



2loi', 



- UJ 2 ) 



and F„, = 



2u 



(w|//-w 2 )' 



(A8) 



where t = 3/(S 3 u±). Here p = (p , p 5 ), p = p aa , 
Pi = Pbb, P2 = Pec, P3 = Re[p 9Q ], Pa = Im^ga], p 5 = F, 
and f = (/o, . . . , f§). The components of f are given by 





-P0 + ^P4, 


(B2) 


fl = 


— Gb Pi + G a b Po, 


(B3) 


h = 


—G c P2 + G ac po, 


(B4) 


h = 


1 

~-p 3 - A/54, 


(B5) 


h = 


~\pi + A P3 ~ t|(2po+Pi+P2-1), 


(B6) 


h = 


+G b pi, 


(B7) 



where and A are given by Eqs. (31) and (32), respec 



tively. To derive Eq. (Bl ), we eliminated the population 
p gg of the ground state from the OBE using the con- 
servation of probability: ^2 n p n n — E The fluorescence 
defined in Eq. ( 34 ) is given by 



Appendix B: System of equations used in numerical 
model 



For computational purposes, 
Eqs. Q and ([34} as 



it is convenient to express 



dp 



f(P,7"), P(-T ) = 0, 



(Bl) 



F(tu)=p 5 (T ;tu). (B8) 
Thus the fluorescence can be obtained by numerically 



solving the system of equations (Bl ), as described in the 
text. 
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